Astronomy & Astrophysics manuscript no. 080923 paper hal 


© ESO 2008 


September 23, 2008 





Extrasolar planets and brown dwarfs around A-F type stars ^ ^ 

V. A planetary system found with HARPS around the F6IV-V star HD 60532. 

M. Desort^ A.-M. Lagrange^ F. Galland^ H. Beust^ S. Udry^, M. Mayor^, and G. Lo Curto^ 

' Laboratoire d'Astrophysique de Grenoble, UMR5571 CNRS, Universite Joseph Fourier, BP 53, 38041 Grenoble Cedex 9, France 
~ Observatoire de Geneve, Universite de Geneve, 51 Chemin des Maillettes, 1290 Sauvemy, Switzerland 
' European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile 

Received date / Accepted date 

ABSTRACT 

Aims. In the frame of the search for extrasolar planets and brown dwarfs around early-type stars, we present the results obtained for 
the F-type main- sequence star HD 60532 (F6V) with HARPS. 

Methods. Using 147 spectra obtained with HARPS at La Silla on a time baseline of two years, we study the radial velocities of this 
star. 

Results. HD 60532 radial velocities are periodically variable, and the variations have a Keplerian origin. This star is surrounded by 
a planetary system of two planets with minimum masses of 1 and 2.5 Mjup and orbital separations of 0.76 and 1.58 AU respectively. 
We also detect high-frequency, low-amplitude (10ms"' peak-to-peak) pulsations. Dynamical studies of the system point toward a 
possible 3:1 mean-motion resonance which should be confirmed within the next decade. 
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1. Introduction 

Radial-velocity (RV) surveys have lead to the detection of nearly 
300 planet^during the past decade. They mainly focus on solar 
and late-type main-sequence (MS) stars (^ F7) which exhibit nu- 
merous lines with low rotational broadening. It is often thought 
that planets around more massive MS stars are not accessible 
to radial-velocity techniques, as they present a small number of 
stellar lines, usually broadened and blended by stellar rotation. 
However, we recently showed (IGalland et al. 1 2005 a. Paper I) that 
with a new radial-velocity measurement method that we devel- 
oped, it is possible to detect planets even around early A-type 
MS stars with large rotational velocities (typically 100km s"'). 
Finding planets around such massive MS stars is of importance, 
as this allows to test planetary formation and evolution processes 
around a larger variety of stars, in terms of stellar mass and time 
scales of evolution processes. This approach is complementary 
with the one which intends to detect planets around evolved 
intermediate-mass stars (e.g.. ISato et al.l 2005. Lovis & Mayor 
2007). In this case, close-in planets have been wiped out but the 
stellar variability is in principle less intense. 

We performed a radial-velocity survey dedicated to the 
search for extrasolar planets and brown dwarfs around a 
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volume-limited sample of A-F main-sequence stars with the 
HARPS spectrograph ( Mayor et aT] 2003) installed on the 3.6- 
m ESO telescope at La Silla Observatory (Chile). We monitored 
a sample of 185 MS stars with B -V ranging between -0.1 and 
0.6. From the measured RV jitter, we computed the minimum 
detectable masses with HARPS, and showed that in 100 cases, 
planets with periods smaller than 100 days can be detected, even 
around stars with early spectral types (down to ~0.1 Mjup at 100 
days around slow-rotating late-F stars). Given the data at hand, 
we also provided the achieved detection limits on the individual 
targets ( [Lagrange et al.| 2008). 

In the course of this survey, we identified a few stars which 
RV variations could be attributed to planets. Most of these stars 
are still under follow-up. We present here the detection of a plan- 
etary system around one of these stars, HD 60532. Section 2 pro- 
vides the stellar properties of this star, the measurement of the 
radial velocities and their relevance; we also present a Kepleiian 
solution associated to the presence of two planets. 



2. Stellar characteristics and measurements 

2. 1 . Stellar properties 

HD 60532 (HIP 36795, HR2906) is located at 25.7 pc from the 
Sun (ESA 1997). Stellar parameters such as mass, age, metallic- 
ity, rotational velocity and effective temperature are taken from 
INordstrom et al.l (2004). and the gravity is taken fi'om |Gray et al. 
(2006). Those values are reported in Table [T] They are in agree- 
ment with a spectral type F6IV-V (F6IV in the Bright Star 
Catalo gue, iHoffleit et al.l (I99I). F6V in the HIPPARCOS cata- 
logue ffiSSl 1997). 
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Table 1. HD 60532 stellar properties. Photometric and astro- 
metric data are extracted from the HIPPARCOS catalogue (lESAI 
1997); spectroscopic data are from iNordstrom et al., (2004'). 



Lomb-Scargle periodogram of the velocities 



Parameter 




HD 60532 


Spectral Type 




F6IV-V 


V sin 7 


[kms-'] 


8 


V 




4.45 


B-V 




0.52 


n 


[mas] 


38.9 ±0.7 


Distance 


[pc] 


25.7 


Mv 




2.40 






-4.94 


[Fe/H] 




-0.42 


Tcff 


[K] 


6095 


log^ 


[cm s"^] 


-3.83 


Ml 


[Mq] 


1 44+0.03 


Age 


[Gyr] 


2.7 ±0.1 



2.2. Radial-velocity measurements 

Since February 2006, 147 high signal-to-noise ratio (S /N) spec- 
tra have been acquired with HARPS, with a S/N equal to 310 
on average. Each spectrum is formed by 72 spectral orders cov- 
ering the spectral window [3800 A, 6900 A], with a resolution 
R ~ 115 000. 

The radial velocities and associated uncertainties have 
been measured with a dedicated tool (SAFIR) which uses the 
Fourier interspectrum method described in lChellil (2000) and in 
iGalland et al.i (2005a). The uncertainty is equal to 0.9ms ' on 
average, consistent with the value obtained from our simulations 
(see Paper I), and includes photon noise, instrumental effects and 
guiding errors (fixed at a upper limit of 0.5 m s"'). 

Note: given the star spectral type, we could also derive the 
RV from gaussian adjustements to the CCFs (Cross-Correlation 
Functions), using a mask with a G2 spectral type. The obtained 
values are comptatible within the error bars with the ones mea- 
sured by SAFIR. 

The amplitude of the radial-velocity variations (more than 
120ms"') as well as their standard deviation (ctrv = 27ms"') 
are much larger than the uncertainties. We now show that these 
variations are not due to line-profile variations and we interpret 
them in terms of the presence of planetary companions. 

The periodogram (Fig.[T]i of the RVs indicates four main pe- 
riodicities in the data. Two of them (near 30 days) are aliases of 
data (Fig. [T] bottom represents the alias generated by our sam- 
pling frequency), whereas the other two (~220 and ~500-1000 
days) are real periodicities of the RVs. This points toward the 
existence of two companions in orbit with periods close to those 
one. 



2.3. A Keplerian origin to the main variations 

In a (RV; bisector velocity span) diagram (Fig. |2|i, the bisector 
velocity spans are spread mainly horizontally and are not cor- 
related with the RVs. This argues in favor of a planetary ori- 
gin to the variations rather than stellar activity (spots, pulsa- 
tions). A spot origin can be rejected from several grounds: 1) 
given the star v sin /, much higher than the instrumental resolu- 
tion (8kms"' compared to 3kms"'), cool spots on the surface 
of the star would induce a correlation between bisector veloc- 
ity spans and radial velocities dPesort et al.l 2007). and the RV 
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Fig. 1. Periodogram of the RVs (top). Four main peaks are vis- 
ible. The two near 30 days are observation alias (bottom, peri- 
odogram of the observation epochs), the other two correspond 
roughly to periods of ~220 and ~500-1000 days. These periods 
are not well constrained because of the limited phase coverage 
(our data span only over ~900 days. The middle panel shows the 
periodogram of the residuals to the final keplerian solution. 



variations would have a period similar to the one of the star's 
rotation periocj^ whereas they are actually much longer; 2) also, 
spots able to produce such amplitude of RV variations would in- 
duce detectable photometric variations, whereas the photometry 
given by HIPPARCOS (ESA 1997) is constant with a scatter of 
only 0.004 mag; 3) the star has a low level of activity (no emis- 
sion in the Ca ii K line. Fig. |3] and low log ^^^)j^ 

2.4. Orbital parameters 

The observed RV can clearly not be fitted by one companion. 
The orbital parameters derived from the best two-companion 
Keplerian solution (Fig. |4|i are given in Table [2] The residuals 
dispersion is 4ms"'. A planetary system of two Jupiter-mass 
planets on long-period orbits fits the data very well. Assuming 
a primary mass of 1.44 ]VI0, this leads to a system with a 2.5- 
Mjup planet on a 604-day orbit and a 1.0-lVljup planet on a 201- 
day orbit. 



- Here, the rotational period of the star P^^t is smaller than ~7 days. 
^ See online section for a figure showing the bisectors of all the spec- 
tra and for an example of an active, not young F6V star. 
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Fig, 2, Bisector velocity spans versus RVs. The bisector velocity 
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Fig. 3. No emission in the Ca ii K line for the HD 60532 spectra. 
This is the average spectrum of all the spectra used after recen- 
tering by the RV variations. 

Table 2. Best orbital solution for HD 60532. 
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Planet b 


Planet c 
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[days] 


201.3 ±0.6 


604 ±9 


To 


[JD-2450000] 


3987 ± 2 


3723 ± 158 


e 




0.28 ± 0.03 


0.02 ± 0.02 


0} 


[deg] 


-8.1 ±4.9 


-209 ± 92 
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[m s-'] 


29.3 ± 1.4 


46.4 ± 1.7 
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^^meas 
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[10-' Mq] 


0.46 


6.2 


Ml 
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Fig. 4. Top: HARPS radial velocities and orbital solution for 
HD 60532, the panel below shows the residuals to the fitted 
orbital solution. Middle and Bottom: Phased fits on individual 
planets. 



2.5. Interpretation oftiie residuals 

The residuals to the possible orbital solutions show variations, 
with standard deviation of 4.4 ms-'; we do not find any peri- 
odicity in these variations. The amplitude of variations is the 
same as the level of the low-amplitude variations seen in the 
line profile. Also, in a (O - C; O - C bisector velocity span) 
diagram, the O - C bisector velocity spans are not horizontally 
spread, which shows that these remaining variations are not due 
to the presence of a smaller-mass planet (Fig. |5]l. They are then 
probably due to stellar intrinsic phenomena. We monitored the 



star continuously for one hour (28 spectra) and measured RV 
variations with a peak-to-peak amplitude of 9 ms-'. This cor- 
responds to the observed semi-amplitude of the residuals. Such 
high-frequency variations argue for stellar pulsations and if we 
would like to increase the measurement accuracy we would need 
to average this eff'ect. Finally, we could not find any periodicity 
associated with the residuals. 
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Fig. 5. (O - C; O - C bisector velocity span) diagram: the O -C 
bisector velocity spans are not horizontally spread (Sect. 2.5 i. 



3. Dynamical issues 

The two orbital periods derived from the fit suggest a possible 
3:1 mean-motion resonance between the two planets. Here we 
investigate the dynamical behaviour of the system to possibly 
distinguish between resonant and non-resonant configurations. 
As we show below, the uncertainty on the semi-major axis values 
for planets b and c (mainly on Planet c, see Table [2]i does not 
permit to state whether the planets are actually locked in mean- 
motion resonance of just next to it. 

The two planets are said to be in mean-motion resonance 
when their orbital periods achieve a simple rational ratio. More 
generally, they will be assumed to be locked ina p + q : p, where 
p and q are integers. \q\ is called the order of the resonance, and 
denotes the number of stable conjunction positions. In the 3:1 
case, we have q - 2. The dynamics within a resonance is usally 
described via the use of the following variable called critical 
argument of the resonance: 



cr - Ac Ab- 



(1) 



where Ab and Ac are the mean longitudes of the inner and the 
outer planet, respectively, and where vJb is the longitude of peri- 
astron of the inner planet (here the less massive one). Resonant 
orbits are characterized by a libration of cr around an equilibrium 
position cTg, while non-resonant orbit exhibit a circulation of cr. 
In non-resonant configurations, the semi-major axes are nearly 
secular invariants, due to independent phase averaging over both 
orbits, while in resonant configurations, the cr libration induces 
secular eccentricity and semi-major axis oscillations. 

We integrate the 3-body system (star + 2 planets, assuming 
sin i - \) starting from the today fitted orbital solution, using the 
symplectic N-body code SyMBA dPuncan et al . 1998). The use 
of a powerful symplectic code allows to adopt a large timestep 
while remaining accurate. We adopt here a timestep of 0.025 yr, 
i.e., a bit less that 1/20 of the smallest orbital period. This is a 
standard prescription in the use of this code to ensure a relevant 
accuracy. We also made tests using a ten times smaller timestep 
to check the validity of our integrations. From a dynamical point 
of view the orbital solution as given in Table |2] is not complete. 
The longitude of nodes Q.b and Q.c are not constrained. Note that 
thanks to rotational invariance, only the difference - ilj, is a 
relevant parameter It will be treated as a free parameter in our 
integrations. 

We first focus on resonant configurations. For this we assume 
for the semi-major axes ab - 0.759 AU and a^. = 1 .57928 AU. 
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Fig. 7. Evolution of the 3:1 critical argument cr over 1000 yr, 
in the same condition as described in Figure |6] We note the cr- 
libration characteristic for resonant motion. 



This choice for a^. (within the error bar) ensures an exact reso- 
nant configuration. 

Figure |6] shows the orbital evolution over lOOyr of the two 
planets in this configuration under their mutual perturbations, for 
a choice of Qf-Q/, = 50°. We note that both the semi-major axes 
and the eccentricities exhibit a superimposition of two types of 
variations: we first have a long period oscillation with ~ 70 yr 
period, in addition to high frequency (1-2 yr periodicity), but 
much smaller amplitude changes. The former oscillation is re- 
lated to the resonant libration motion (see Fig.[7]i, while the lat- 
ter is phased with the synodic motion of the two planets. It is 
thus related to the mutual perturbations of the two bodies at con- 
junction. The fact that this high-frequency term shows up here 
is an indication for weak - but significant - chaos in the system. 
The integration was extended up to 10** yr. Over this time span, 
the planetary system does not show any indication for instabil- 
ity. The semi-major axes of the two planets oscillate between 
0.748 AU and 0.767 AU for Planet b, and between 1.558 AU and 
1.6 AU for Planet c. The eccentricity of Planet b oscillates be- 
tween 0.13 and 0.32, and that of Planet c between and 0.14. 
Note that the present value of Planet c's eccentricity is close to 
the bottom of its actual variation range. 

Figure [6] shows that, given the error bars listed in Table [2] 
the secular evolution of the semi-major axis of Planet b should 
be detectable within ~ lOyr from now. In fact, what is expected 
to be detected is only the low frequency secular variation. The 
high frequency oscillation remains below the error bar 

Now, we point out that Fig.|6]corresponds to Q.c-Q.b - 50°. 
Figure [8] shows the evolution of the semi-major axes in the same 
conditions, but starting with Qc - Q/, = 150°. We still see the 
same kind of oscillations, but phased differently. In Fig. [6] the 
semi-major axis of Planet b increases in the first 10 years, while 
in Fig. [8] it decreases. 

Although the actual value of Q^. - Q.b is unknown, the de- 
tection of semi-major axis evolution of Planet b within 10 years 
should be considered as a strong indication of a resonant config- 
uration. Figure |9] shows the evolution of the semi-major axis of 
Planet b and of the critical angle cr for 3:1 resonance for an ini- 
tial choice of flc = 1 -57528 AU. As can be seen from the evolu- 
tion (circulation) of cr, this configuration is non-resonant, while 
still within the error bar of the orbital fit. Hence from the or- 
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Fig. 6. Orbital evolution over 100 yr of the semi-major axes (top) and eccentricities (bottom) for planets b (left) and c (right), under 
their mutual perturbations, in a 3:1 resonance configuration, for an initial choice of Q.^ - Q./, - 50°. 
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Fig. 8. Evolution of the semi-major axes of the two planets in the same condition as described in Fig.|6j except that we assumed 
Q,. - Qfe = 150°. 
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Fig. 9. Evolution of the semi-major axis of Planet b and of the critical angle cr for 3:1 resonance in the same conditions as described 
in Fig.[6j but with Oc - 1.57528 AU. This is a non-resonant configuration. 



bital fit itself, it is not possible to definitely state whether the 
two planets are actually locked in 3:1 resonance or not. From 
Fig. [9] we see that in non-resonant configuration, the variations 
of the semi-major axis of Planet b achieve a much lower ampli- 
tude than in the resonant case. This remains true for any choice 
of initial non-resonant configuration. In fact, in that case, we 
only have the high-frequency, small-amplitude oscillation. As 
expected in non-resonant case, there is no secular evolution of 
the semi-major axis, contrary to the resonant case. The net con- 
sequence of this is that given the error bars of the fit, no semi- 
major axis variation for Planet b should be detected within 10 
years in the non-resonant case. 

Our conclusions concerning this dynamical study are thus 
the following: 

- Given the error bars of the fit, it is not possible to definitely 
state the planets are actually locked in mean-motion resonant 
or not; 

- The two-planet system is dynamically stable although sig- 
nificantly chaotic; 

- Although it is impossible to predict the variation sense due to 
unconstrained 0.^-0.1,, any detection of variation in the semi- 
major axis of Planet b within 10 years from now should be a 
strong indication for a resonant configuration. Further mon- 
itoring of this system should therefore be initiated to detect 
this possible variation. 



4. Concluding remarks 

We have shown that HD 60532, an F6IV-V, 1.44Mo star hosts 
two planets with minimum masses of 1 and 2.5 Mjup and orbital 
separations of 0.76 and 1.58 AU respectively, in a possible 3:1 
resonance which needs to be confirmed within the next 10 years. 
Noticeably sofar, only one other multiple system had been re- 
ported around a MS star more massive than 1.3 Mq (HD 169830; 
1.4 M0). Note also that the low metallicity of HD 60532 is not 
common for stars harbouring Jupiter-mass planets; the relation 
between the star's metallicity and the presence of massive plan- 
ets has still to be investigated further 
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Appendix A: Additional material 



Figure [AT| presents the bisectors of the CCFs (which stands for 
all the lines) for the spectra from which we measured the radial 
velocities of HD 60532. It shows that spectra are mainly shifted, 
with yet small line-profile variations superimposed to the shifts, 
and thus that the radial velocities measured are not induced by 
line-profile variations. 

Figure A. 2 shows an example of an active F6V star 
(HD 30652, V sin / - 16kms ') that we observed during the same 
survey, and for which RV measurements are correlated with line- 
profile variations, indicating that they are are induced by a spot. 



List of Objects 

'HD 60532' on page[Tl 
'HD 169830' on page]6] 
'HD 30652' on page|6f 
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Fig.A.l. Bisectors of the CCFs. 
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Fig. A.2. Example of an active F6V star for which measured RVs 
(top) which are correlated to the bisector velocity spans {bot- 
tom). 



